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We study current correlation functions in a diffusive junction out of equilibrium. We calculate 
corrections to the electric current and to the zero frequency shot noise due to electron-electron 
interactions. Contrary to the equilibrium situation (where the corrections to the current and to the 
current noise are related through the fluctuation-dissipation theorem (FDT)), these two quantities 
behave differently: the correction to the electron current are governed by the largest of temperature 
and applied voltage, while the correction to the shot noise is always governed by the temperature. 
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I. INTRODUCTION 

The physics of non-equilibrium mesoscopic systems 
was thfij^ubject ofj^n extensive research for the last 
decadeLfd. KhlusQ, who employing a semiclassical 
approximation considered a multichannel two-terminal 
junction, was the first to notice that the non-equilibrium 
current noise in such a system strongly depends on the 
transparency of the barrier, and vanishes in the limit of 
completely transparent junction. This is in contrast to 
the conductance which remains finite in this limit. Later 
LesovikQ has shown that for coherent transport through 
a multi-channel two-terminal system the zero frequency 
current noiscO is given by 

S(0) = ^ / dEj2{TniE)[fL{E)il-ME)] + 

fnm^~fL{E)]-T^{E)[fUE)~fR{E)r} . (1) 

This was derived using the Landaueifl scattering states 
approach. Here / is the grand-canonical Fermi-Dirac dis- 
tribution function, L and R denote the left and right 
reservoirs respectively, taken at different values of the 
chemical potential, and {Tn} are the channel transpaijen- 
cies. This result was later generalized by BiittikerO to 
include the multi-channel multi-terminal case and exper- 
imentally confirmed by Reznikov et. alB. The zero tem- 
perature limit of eq.(|I]) yields the shot noise: 

^i^)^i^J2^"{l-Tn)eV . (2) 

n 

The remarkable dependence of the shot noise on the 
quasiparticde. charge ( as was first emphasized by Kane 
and Fishediy) allowed the first direct measurement of 
fractiopal charge in the context of EDHE by Reznikov 
et. al]13 and by Saminadaya et. alx3 The suppression 
of the shot noise for open channels ( in full accordance 
with the semiclassical resulto, eq. (||)) is a manifestation 
of correlations among the electrons, arising from Pauli 
exclusion principle. In disordered systems one usually 
cannot control (or obtain information about) the indi- 
vidual channels; instead it is useful to consider average 



quantities. For coherent diffusive systems (smaller than 
the inelastic, dephasing and the localization lengths) the 
various moments of the current-current correlation func- 
tion have boeii calculated employing random matrix the- 
ory (RMT)t3. Interestingly enough, it turned out that 
the fermionic suppression of the noise present in ballistic 
systems is manifest in the configuration averaged noise 
(in diffusive systems) as well: the shot noise is suppressed 
by a factor of 1 /3 as compared with the result expected 
for classical particles. 

Other than the Landauer scattering states approach, 
one may use the kinetic equation, or the so-called Kogan- 
Shulman approachEj. It is valid for systems where the 
dynamics (but not necessarily the statistics) of the parti- 
cles is classical. Comparing with the scattering states 
approach, the latter is not restricted solely to non- 
interacting electrons. Unfortunately, it is limited by the 
usual conditions of applicability of the kinetic equation: 
it does not yield quantum interference effects, and can- 
not be used for calculation of the cummulants higher than 
two. To overcome the lattejt-shortcomings one may em- 
ploy the Keldysh formalismll3. 

In the present analysis wc find it convenient to employ 
this formalism in the form of a non-linear sigma-model. 
This machinery was found to be an effective tool for de- 
scribinft-the low energy physics of disordered electronic 
systemsEEl. We will not attempt to review here the various 
technical steps involved in the derivation of the model, 
but rather recall some of the main ingredients involved. 
So long as non-interacting electrons are concerned, the 
super-symmetric sigma-modelll3 turns out to be a par- 
ticularly suitable tool for calculation. The inclusion of 
electron-electron interaction appears, though, to render 
this approach infeasible. Interaction among the electrons 
could be included within the the replica sigma-model, 
which has been introduced by Finkelstein for the study of 
disordered Fermi liquidsta. As haSjXecently been demon- 
strated by Kamenev and AndrjeevEfl and independently 
by Chamon, Ludwig and Nayakc^, there is an alternative 
approach of constructing a sigma-model for interacting 
electrons - the so-called Keldysh sigma-model. 

To wrap up this introduction, we find it useful to de- 
fine some of the physical regimes to be discussed below. 
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We first stress that the notion of high/low frequency car- 
ries two meanings. As far as the diffusive motion of the 
electrons (viewed as independent particles) is considered, 
one can refer to the frequency as "low" if it is smaller 
than the Thouless energy of the system, Eth ~ hD/L^. 
As we study the case of interacting electrons, another 
criterion arises, relating to the propagation of electro- 
magnetic fields in the system. Consider the simplest sit- 
uation when the effects of external screening (by, e.g., ex- 
ternal gates) can be neglected. It is readily seen that the 
solution of the Maxwell equations (or the self-consistent 
solution of the Maxwell equations combined with the ki- 
ef. ~ 



III below) depends on the 
For a standard 



netic equations, cf. Section 
effective dimensionality of the sample 
three dimensional geometry the frequency at which the 
current noise is measured (w) should be compared with 
the inverse Maxwell time (ujm — Ana). At {uj <^ oJu) the 
"cross-coordinate" current correlation (measured at dif- 
ferent cross-sections) is independent of the spatial coordi- 
nates (and their separation, lx,x')- Having this electrody- 
namics in mind, it is this limit which should be referred to 
as "low frequency " . For quasi- two dimensional films this 
regime is realized under the condition: lo <C H2D/lx.x' 
where the inverse screening length K2 = 2'r:e^v. The ef- 
fects of external |Sereening may modify the effective elec- 
tron interaction,EirE3. Our point, though, is that quite 
generally there is some frequency, below which the cor- 
relation function docs not depend on the spatial coordi- 
nates, i.e. it is a constant as function of In the 
case of interacting electrons, in addition to the long wave- 
length propagation of the electro magnetic field, one also 
needs to account for inelastic collisions among the elec- 
trons. We thus introduce another dimensionless param- 
eter, namely the ratio between the system's length, L, 
and the inelastic length lin. 

The outline of the present work is the following: In 
the next section we consider non-interacting spinless elec- 
trons in a weakly disordered conductor neglecting weak 
(and in two-dimensions strong as well) localization ef- 
fects. The latter can be justified once the quantum co- 
herence of the electrons is destroyed by certain dephasing 
mechanisms (such as fluctuating electro-magnetic fields) , 
or is restricted by the finite size of the system. In other 
words, we consider a situation where the localization 
length, ^, exceeds either the dephasing length Iff, or the 
system size (^ ^ maxjZ^jL}). We extend the Keldysh 
sigma-model formalism to the case of systems out of equi- 
librium (including the effect of restricted geometries). In 
Section [II we review the calculation of the current noise 



(in and out of equilibrium), discussing both the low and 
the high frequency limits, paying special attention to 
the non-homogeneity (in space) of the current fluctua- 
tion in the later case. In Section |^ we consider the 
effect of electron-electron interactions. We review results 
obtained from the kinetic equation approach, referring 
briefly to the limits of Im «C L and Im ^ L. We then fo- 
cus on the latter limit presenting results for the effect of 
electron-electron interactions on the noise which go be- 



yond the kinetic equation (eqs. ^ and pq ). Our res 
are analogous to the Altshuler-Aronov correctionsoc3 
found earlier for equilibrium noise. This is done using a 
model short-ranged interaction. The general expressions 
obtained for the noise are then analyzed in the interest- 
ing case of two-dimensional geometries. We show that 
electron-electron interaction leads to the suppression of 
the non-equilibrium current noise which is stronger than 
the suppression of the current itself - a manifestation of 
many-body correlation effects. Eqs. (^, [T^ and |6^) are 
the main results of our analysis. The three appendices 
include some technical parts of the derivation. The anal- 
ysis p«3sented here is an extended account of an earlier 
workEa. 



II. KELDYSH FORMALISM: 
NON-INTERACTING ELECTRONS IN 
DIFFUSIVE SYSTEM 



The theory for quantum systems out of equilibrium was 
suggested] simultaneously by Kadanoff and Baym and by 
KeldyshcZl. This technique has beea widely used on the 
level of diagrammatic calculation£j . Such techniques 
can alsQ-be developed for the study of weakly disordered 
sample^J. Instead of using here a diagrammatic tech- 
nique, we will present a path-integral description of the 
problem. For free fermions the action can be written as 



S= dr dtC , 



(3) 



where the Lagrangian density is given by 

C = ¥[Go ' - C/*s]* . (4) 

Here the Schrodinger operator of an electron of mass m 
in the presence of a vector potential ea/c is given by 



(5) 



The energy is measured from the chemical potential fi. 
The action now contains an integral over space and an 
integral in time over a Keldysh contour (see Fig |l|) . In 
analogy with ordinary diagrammatic techniques one de- 
flnes an electron Green's function, i.e. a time ordered 
correlation function 

G„^(r, t; r\t') = -i{T,^(r, t,,)^\r' ,t'f,)) (6) 

Here is a time ordering operator along a Keldysh 
contour, the indices a, fi denote the branches of the 
Keldysh contour (minus or plus for the upper and lower 
branch of the contour respectively). Here () denotes av- 
eraging with respect to the action eq. (||) . 
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FIG. 1. A Keldysh contour 
It is also useful to introduce the generating functional 



Z[Sl] = / V^exp{iS[^,a]/h) ^ 



(7) 



One notices that the generating functional (|^) is identi- 
cally equal to unity, so long as the external source is the 
same on the upper and lower branchestJ {a-{t) = a+(t)). 
For this reason one may explicitly perform the averaging 
of the generating functional over disorder. Below we con- 
sider short-range and (5-correlated disorder, which leads 
to following definition of disorder averaging: 



1 2 



exp < — [AttUvt) dr 



J VUdis exp ^-TTt/r j dr Uj^^{r) 



(8) 



where r and v are the elastic mean free time and the 
single-particle density of states respectively. When ap- 
plied to the generating functional it produces an effective 
non-linear action 



J dr I dt*G;7V(3^*- 



dr 



where 



1 
1 



(9) 



(10) 



The non-linear term in eq.(|^) can be decoupled with an 
auxiliary bosonic (a 2x2 matrix) field Q{r,t,t') = Qr,t,t', 



VQ exp 



dvdtdt' 



dt^{Y,t)(j^^^-^{v,t) 



At 



2t 



(11) 



As a result of this Hubbard-Stratonovich transforma- 
tion one obtains the following action: 

^S[Q] = -^Tr{Qn + Trln[(^G-i + ^Q^ a^^)] , 

(12) 

It follows from equation (||) that the components of the 
matrix Green function are linearly related 



G— + + G+— 



(13) 



It is thus convenient to perform a rotation in the Keldysh 
space 



Here the unitary rotation matrix R is given by 
1 / . ^9"i\ 1 / 1 — 1 



R = 



V2 

yielding (cf. eq. ( [12| 



V2\^ i; ' 



(14) 



(15) 



iS[Q] 



Tr{g^} + Trln[Go-^ + -Q] . (16) 



Clearly, the resulting action is exactly equivalent to the 
original one, eq. (^). However, it turns out that in many 
physical problems the fast variations of the matrix Q in 



space and time are unimportant, and we may integrate 
them out. In the weak disorder case {epT 3> ti, where 
is a Fermi energy) the main contribution comes from 
the minimum of the action, justifying the use of a saddle- 
point approximation. Varying the action with respect to 
Q, one arrives at the following saddle-point equation 



Go V(o) + 



(17) 



Taking the Fourier transform with respect to t — t' one 
can show that the matrix 



A(r,6) 



1 2F(e) 
-1 



(18) 



( for an arbitrary F) is a solution of equation (g_7|) . Apart 
from F(e) which needs to be determined, the solution 
for the saddle-point (eq. |lj) is not unique. Indeed, the 
action ( p^ ) is invariant under any homogeneous unitary 
transformation 



Q = C/UC/ 



(19) 



To calculate the fluctuations around the saddle 
point (^ we will distinguish between longitudinal and 
transversal fluctuations, the latter do not modify the 
value of Q^. The spectrum of the longitudinal fluctua- 
tions has a gap; for this reason they are irrelevant for the 
low energy physical properties of the system (although 
they do affect the bare values of the parameters). The 
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gapless transversal fluctuations (soft modes) should be 
taken into account. 

We now focus on the long wavelength fluctuations 
around the saddle point A. To study such fluctuations 
we now extend our discussion and assume the matrix U 
in eq. (|l^) to be a slow function of the spatial coordi- 
nate. One may now expand the action in gradients of 
U{r). We thus obtaiu-a new action in a form which is 
usually referred to aa^ a non-linear sigma model 

iS[Q] = -^Tt[D {VQf + AieQ] , (20) 

where the integration over the field Q is now restricted 
by the non-linear constraint 



= 1 . 



(21) 



One may detive an equation of motion ( quantum ki- 
netic equatiorJla) , which is compatible with eq. (^): 

DV{QVQ)+i[e.,Q]^0 . (22) 

We now use the following parameterization: 



Q = Aexp(W^) , 



(23) 



where the matrix Wx^e,e', in turn, is parameterized as 
follows 

\ —Wx,e,e' —Wx,e,e'^x,e' J 

and where the fields -F, w and w are unconstrained. Note 
that A is parameterized by F as well, (cf. eg. ([l8|) ) . Ex- 
panding the action in soft mode fluctuaticms around the 
saddle point we obtain, to second order inE2l w,w: 

iS[W] = iSo[W] + iSi[W] + iS2[W] . (25) 

The zeroth order term is identically equal to zero 

iSo[W] = 0. (26) 

This implies that, in the absence of external sources, the 
generating functional is equal to unity even on the level 
of the saddle point equation. The term linear in W is 
given by 

iSi [W] = TrvhTT{ie{wF - Fw) + DwV^F}. (27) 



Selecting A to be a saddle point of the action ( |20| ) , which 
requires the part linear in W to vanish as well, we obtain 



^{el-e2)Fx,e^,c,- DV^Fx,e^. 







(28) 



The function F{x,e) is related to the single particle dis- 
tribution function /(x, e) 



Indeed, eq. (pSD coincides with the Boltzmann equation 
within the diffusion approximation for /. However, since 
eq.(|2^) is a linear homogeneous equation, F may be re- 
lated to / only up to multiplicative factors. To establish 
the relation, eq.(|29|), one needs to resort to the inhomo- 
geneous version of eq (^), where the r.h.s. is replaced 
by a collision integral due to electron-electron interac- 
tions. This indeed has been done in Ref.E3 ( cf. with 
their eq.(138)). 

We can describe what we have done above slightly dif- 
ferently. We note that the parameterization, eq.( p3|), i s 
a particular choic^ satisfying the constraint, eq.( 
We may look for a particular case for which A (which, 
to begin with, was a saddle point of the original action, 
eq. ([iq)) is a saddle point of the non- linear cr- model ac- 
tion, eq. ( |20| ) (i.e., S'i[W^] = 0). This leads to equation ( 
(p8|)). On the other hand , we know that one can employ 
the Keldysh technique perturbatively, and relate A (with 
the parameterization of eq. (|23|)) to G(r, r) ( where G 
is the Green's function in the Keldysh space), with the 
(1, 2) element thereof related to the distribution function 
/. This justifies referring to eq.(p8|) as a kinetic equa- 
tion. To consider larger deviations from the saddle point 
i^one needs to consult the exact saddle point equation 

Equation (|2^) (or ^8|) should be supplemented by ap- 
propriate boundary conditions at the endpoints of the 
sample. The quadratic fluctuations around the saddle 
point are given by 



Wx,<i,€'[~DV'^ + i(e - e')]wx,e',e 



(30) 



This action describes the dynamics of the fluctuating 
fields w and w. The dynamics may be affected by the 
last term on the r.h.s. of eq.(^o|), accounting for devia- 
tions from thermal equilibrium. 

We now specialize to a concrete geometry and specific 
boundary conditions. We consider a system made of two 
clean ( i.e., no disorder) and stationary reservoirs, con- 
nected to each other through a diffusive bridge as shown 
in Fig. (I). 



^l+eV/2 



-eV/2 



A 



F{x,e) = l-2S{x,e) 



(29) 



FIG. 2. A diffusive bridge (the shaded area) between 
two reservoirs 
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non-zero bias {eV ^ 0) gives rise to a finite d.c. current 
through the bridge. Hereafter we assume that the typical 
energy relaxation time is much larger than the time-of- 
flight of the electrons through the micro-bridge, resulting 
in the following stationary diffusion equation 



DW'^F{r,e) = 0. 



(31) 



Note that under the conditions specified above, eq. ( |3l| ) 
contains no energy collision integral and no time deriva- 
tive. Furthermore, we assume that the typical time a 
particle spends in the reservoir is much larger than the 
energy relaxation time, so that the single-electron distri- 
bution function is subject to local equilibrium boundary 
conditions at the edges, with the respective chemical po- 
tentials as indicated in Fig. (||). It turns out that the 
solution to. eq. consistent with such boundary con- 
ditions ig2il 



F(e,x) = (|)Fo(e+ey/2) + 



l--jFo(e-ey/2). (32) 



Here 



Fo(e) - tanh (^) 



(33) 



As we have prepared the tools for our analysis of non- 
equilibrium noise, we are now set to construct the corre- 
lation functions for w and iJu. It is convenient to define 
the diffusion propagator as 

{-iw + DV^) D(r, r', uj) = ;^'5(r - r') , (34) 



with the boundary conditions defined on the edgesOu. 
Since no current may flow through infinite walls , the 
derivative of the diffusion propagator must vanish in the 
transversal direction, dD/dy — dT)/dz — 0. Moreover, as 
the diffusion modes cannot propagate through the clean 
metallic leads, the diffusion propagator must vanish at 
the contacts to the latter. 

The action (|o|) can be viewed as a differential oper- 
ator acting in the 2x2 space spanned by the fields w 
and w. We need to invert this operator to obtain their 
correlation function. 

(w(a;,ei,e2)w(x',e3,e4)) = 

2(27r)2,5(ei - ei)5{e2 ~ e3)D(x, x', ei - £2) , (35a) 
{mix, ei, £2)^(2;', £3, £4)) = -(27r)^(5(£i - ei)5{e2 - £3) • 

(35b) 

(tD(a;,£i,£2)w(a;',£3,£4)) = , (35c) 

were the dimensionless conductance g ~ hvD. Note 
that as a result of the Keldysh rota tion, eq.(|lj), the 
correlation function oi {ww) (eq.(35c)) vanishes. Next 
we employ the action ( |30| ) to calculate various observ- 
ables. For this purpose we will use the generating func- 
tional, eq. (0). To write it in terms of the Q matrices 



we note that the magnetic field enters the action ( |2C 
through the "long derivative" d — V + z[a"7Q,, ], where 
ea/c is the vector potential. This follows from the gauge 
transformationE^ 5* Ua^ which allows to eliminate 
the vector potential from the action. Provided that the 
spectrum of the electrons can be approximated as linear, 
one chooses VJ7a = —iAtUa- As a result the generating 
functional (^ can be written as 



Z[a\ = j'DQ cxp (^-^Tt[D {dQf + ^Q] 



(36) 



As a simple example for an observable we consider the 
mean electron current through the bridge (at a vanishing 
magnetic field), which is given by 



-e5Z[si] 



2i 63.2 



(37) 



a2=0:ai=0 



The corresponding leading diagram for the average cur- 
rent is shown in Fig.(||) 




G 



K 



FIG. 3. The leading diagram for the mean value of 
the electron current in the diffusion approximation. The 
shaded semi-circle represents the Keldysh component of 
the Green function. 



Employing eq. ( p6D we find an expression for the d.c. 
current 



/ = envnDi:T{{QVQ - (VQ) Q)72} 



(38) 



Taking the saddle point solution for Q, eq. (18), one 
obtains 



^Ohm = Gohm^ , 



(39) 



which is Ohm's law; for diffusive systems of length L and 
cross-section A the low frequency conductance is equal 
to 



Gohm ~ ^~[, ' 



(40) 



where the conductivity a is related to the diffusion coef- 
ficient and the thermodynamic density of states through 
the Einstein relation 



(41) 
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III. CURRENT NOISE IN DISORDERED 
JUNCTIONS 

Besides the average value of the current one may study 
its fluctuations in timallEZl. As far as the second mo- 
ment is considered, one is oftpa concerned with the sym- 
metrized correlation functiorH 

S(t, X, t', x') = i((/(t, x)i{t', x') + /(f, x')i{t, x))). (42) 

In a steady state S" is a function of the difference between 
its respective arguments. In equilibrium this quantity 
is related to the linear conductance by the Fluctuation- 
Dissipation Theorem (FDT) 



S''*(w) = ?icjcoth 



flLO 

2T 



G{u), 



(43) 



where G{uj) is the frequency dependent conductance. 
Note that the choice of the microscopic parameters ( 
such as impurity strength and concentration ) enter only 
through the diffusion coefficient, hence the conductance. 
The relation, eq. (^ ) , between the conductance and cur- 
rent noise is universal at equilibrium and does not depend 
on the microscopic details of the system at hand. Out of 
equilibrium there is no fundamental expression govern- 
ing the relation between conductance and the spectral 
function of the noise, thence one needs to calculate them 
separately. 

To calculate a two-operator correlation function, one 
may choose the time indices to relate to different 
branches of the Keldysh contour, and consequently repre- 
sent the correlator as a Keldysh time-ordered expression, 
making it possible to apply Wick's theorem. Formally, 
{A{t)B{t')) = {TcA{t+)B{t'-)), where A,B, are arbi- 
trary operators and Tc denotes the time ordering opera- 
tor on the Keldysh contour. The symmetrized current- 
current correlation function (at a vanishing magnetic 
flux) can be represented as 



S{x,t;x',t') = 



4 6a2ix,t)dsL2{x',t') 



■ (44) 



ai=0;a2=0 



Here we have used the identity Z[ai;a2 — 0] = 1. After 
functional differentiation one obtains the following equa- 
tion 



S{x,t;x',t') = 



\ x.t\x^t 



TTVD 



M,^tM,,,t') , (45) 



M{x,t) ^ Tr<^ / dti 



Q{x,t,ti)VQ{x,ti,t) 



{WQ{x,t,ti))Q{x,ti,t) 



72 



(46) 



One may note the correspondence between the expres- 
sion ( [pt ) and the direct diagrammatic analysis published 
earlieiEa (cf. Fig. 0). 






FIG. 4. The leading diagrams for the current noise; 
the shaded semi-circle represents the single particle dis- 
tribution function; the dashed lines represent contrac- 
tion of impurity lines (i.e., ensemble averaging) 

The first (diamagnetic) term in eq. ( ^ ) is local in 
space (i.e. decaying on distance of the order of the elas- 
tic mean free path) and corresponds to the bare diagram 
(^) , which describes the fluctuations of the distribution 
function at a given point. The non-local part corresponds 
to the diagrams (^a-f), which describe the diffusive prop- 
agation of current fluctuations throughout the sample. It 
is convenient to define a local effective noise temperature 



,iuj,x) 



dell 



F , F ] 



(47) 



such that at equilibrium the noise can be cast as purely 
"thermal" : 



Ceg _ r,rp 



(48) 



One may work out (see Appendix ^ the explicit form of 
the current-current correlation function: 



S{x, x , uj) 



,Xx,u!)S{x — x') + 



where the subscript denotes averaging with the action 
(pO|); we have also introduced the notation 



(49) 



1 
1 



1 

1 



The same expiiession was independently derived by Na- 
gaev recentlja. 



^x.t t' — ^AQx,t.,t'l2Qx' ,t' .tl2 — 5t^t'Qx.,t,t'Qx' ,t' ,t \5x,x' 



Although eq.(49) is generally quite complicated, one 
may study certain limiting cases thereof. 

In the low frequency limit ( uj ^ ETh) the correla- 
tion function does not depend on the spatial coordinates 
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due to particle conservation. To study the temperature 
and frequency dependence it is convenient to integrate 
eq. ( ^Of ) over the volume of the system , which leads tot3 

S{lo) = ^[4S"«(w) + S"'(w + eV) + S'^uj - eV)] . (50) 
6 

In the high frequency limit {uj 3> Exh) the behavior 
is qualitatively different. The "cross-coordinate" noise 
spectrum (cf. eq.(|49|)) in this situation is a decreas- 
ing function of the distance between cross-section {x, x'), 
through which the current fluctuation are studied. If the 
distance Ix.x' = |a; — is much larger than the typical 
diffusion length y/{D/uj), the correlation function ( |49| ) 
practically vanishes. 

We stress that while the electrons here are assumed 
to be non-interacting, the low frequency limit of equa- 
tion ( |50| ) is valid for interacting electrons as well (in the 
limit of low electron-electron collision rate) , and is in 
good agreement with experimentcH. To see this onCpUiay 
follow arguments similar to those given by NagaevEllEj. 
Let us assume that it is possible to describe the kinetics 
of interacting electrons by a Kogan-Shulman type set of 
equations 



6j = -eDVSn + a{oj)5E{uj) 
iueSn + "S/Sj — , 
V5i^ = ineSn . 



(51) 



Here the fluctuations of the electric current density, Sj, 
are related to the fluctuations of the particle density , 
6n, and the electric field, SE, in a self-consistent field ap- 
proximation. For simplicity we have assumed here that 
the elastic relaxation time and the electron density of 
states are constant in the vicinity of the Fermi surface. 
Furthermore, we have ignored the contribution to the re- 
sistance associated with elastic scattering at the contacts 
to the leads. If this is not-the case, the eq.(pTl) should be 
modified correspondingly^ 

The "random source" 6j^^* is generated by short-range 
scattering events of conducting electrons on the impuri- 
ties. Motivated by the fact that scattering events on 
scales larger than the mean free path are uncorrelated, 
we model the random source correlation function by 

{Sf:'Hr)Sjir'ir'))u = ^7:a{u)6ir - r')T{r,oj) . (52) 

In the low frequency limit the continuity equation guar- 
antees that the current-current correlation function is in- 
dependent of the choice of cross-sections, i.e., it is con- 
stant as function of its coordinates. To find this constant 
it is convenient to integrate the current correlation func- 
tion with respect to its coordinates over the volume of 
the sample. In our model there are no fluctuations of the 
electron density and electrostatic potential in the leads. 
This implies that the gradient term disappears after the 
volume integration, i.e., we rederive the non- interacting 



result, eq. (pO[) . Given the fact that neither the diffu- 
sion propagator nor the effective noise temperature en- 
tering eq. ( ^ ) depend on the phase of the electron wave- 
function, we conclude that the current noise in diffusive 
systems is unaffi 
in the systemoE 



iited by quantum decoherence processes 
to leading order in 1/g. 



IV. CURRENT FLUCTUATIONS OF 
INTERACTING ELECTRONS 

So far we have considered the physics of non- 
interacting electrons. In the present section we study 
how the electron-electron interaction affects the current 
through a diffusive junction as well as the current fluctu- 
ations, all this out of equilibrium. One may note that the 
conductance ( including weak localization corrections) is 
insensitive to the shape of the electron distribution func- 
tion. By contrast, the double-step shape of the distribu- 
tion function is crucial as far as shot-noise is concerned. 
This provides us with a motivation to consider the effect 
of electron-electron interaction on the single-electron dis- 
tribution function, hence the noise. Since the electrons 
are interacting particles, one should include the (electron 
density dependent) electric field in the kinetic operator, 
supplemented by a self-consistency condition, eq.(pT|). 
This leads to the appearance of a smooth electric field 
in the constriction which, on its own (as we have already 
seen in the Section ^, does not give rise to any effect 
on the low frequency noise, but may modify the higher 
frequency spectrum of the current noise. This case has 
been studied extensively applying tbe|-Sb«lman-Kogan 
approacht3to various physical systcmEJoEil. Briefly, the 
presence of interaction defines a characteristic RC time 
scale. Depending on the ratio between the measured fre- 
quency and the inverse of this time scale, fluctuations of 
the distribution function propagate (or do not propagate) 
through the entire system ( over time scales l/oj ). In 
the low frequency limit, fluctuations of the distribution 
function do propagate throughout the entire system (over 
time scales ~ l/"-^), in which case current-current corre- 
lators are independent of the spatial coordinates and do 
not discriminate between interacting and non-interacting 
electrons. The value of the shot noise is then given by 



S(0,^f 



(53) 



By contrast, at high frequencies ( larger then the inverse 
RC time) interactions affect thej-frequency and the space 
dependence of the current noiseEj'O. 

However, in addition to the creation of a smooth elec- 
tric field, electron-electron interactions may lead to in- 
elastic scattering. The relative importance of such a 
scattering process is determined by the ratio between the 
sample size L and the inelastic mean free path When 
the inelastic length is much larger than the system size, 
electron-electron scattering leads to a small positive cor- 
rection to the noise given by eq.(pO[) (cf. Rcf.(Ell)). For 
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systems with a short relaxation length ^ L), the 
distribution function can be approximated by a quasi- 
equilibrium one. The latter implies that each piece of 
the conductor possesses its own effective temperature and 
electro-chemical potential. The contributions of the in- 
dividual pieces are then added independently. The local 
quasi-equilibrium conditions lead to the "smearing" of 
the double-step distribution function obtained for non- 
interacting electrons out of equilibrium, eq. , and con- 
sequently to the enhancement of the current noise (or 
in other words to ar^iartial suppression of the electron- 
electron correlation )li2lc3 



S(0) 



el 



(54) 



We now focus on the limit of large inelastic length, 
and study the effect of electron-electron interaction be- 
yond the applicability of the kinetic equation; we evaluate 
the interaction corrections to the current noise. While 
this contribution is small in magnitude, for system with 
effective dimensionality d < 2 it has strong (singular) 
dependence on temperature and voltage, contrary to the 
large ( but practically constant ) non- interacting part. 

We take the function F to be that of the non- 



interacting gas (in general out of equilibrium), eq.(32). 
Corrections to the equilibrium conductance ( and noise) 
due to the interplay between interaction and disor- 
der ha ve ,, o riginally been discussed by Altshuler and 
AronovESnJ and in two-dimensional settings are given by 



SG 



G 
2^g 



In 



Tt 
IT 



(55) 



One may note that unlike the non-interacting case 
(cf. eq. (|4C|)), the conductance due to the Altshuler- 
Aronov term depends strongly on temperature, i.e. on 
the smearing of the distribution function. In the non- 
equilibrium problem the shape of the distribution func- 
tion is markedly different from the Fermi-Dirac distri- 
bution function, which motivates us to study analogous, 
out-of-equilibrium corrections to the current and to the 
noise. In order to do so we incorporate the Coulomb 
interaction into the sigma model. The action of the in- 
teracting electrons may be written as 



(56) 



where 5*0 [^'] is given by eq. (j|) and the interaction part 
of the action can be described as 

S^ntm = -l^ / drdr'dtp.,{r,t)afVo{r-r)p,{r',t). (57) 
1=1-' 

Here p is the electronic density, pi = while Vb rep- 

resents the interaction potential. One may introduce a 
new auxiliary bosonic field 



(58) 



which decouples the interaction in the particle-hole chan- 
nel, yielding 



iS[^ = i?iTr{$'^t/o"V(3)|.} . 

z5[|.,Q]=-^Tr{Qn + 
4r 



Trln 



G,^a 



1^(3)- 



2t 



(59a) 
(59b) 



(59c) 



A major difficulty now arises. The study of an electron 
gas at equilibrium relies on the fact that the correlation 
function of the fluctuations of the electric potential in 
the sample can be expressed through the polarization 
operator of the gas. This is not the case out of equi- 
librium. To find such a correlation function in the non- 
equilibrium case is a demanding task, which, in its gen- 
eral formulation, is yet to be solved. Formally the diffi- 
culty arises because the conventional RPA approximation 
has to be modified. Physically this is related to the fact 
that screening of the interaction for the non-equilibrium 
gas differs from the standard equilibrium problem. We 
circumvent this difficulty by considering a toy model of 
an instantaneous short-range interaction: 



Vo{r-r';t-t') = hrS{r-r')S{t~t') 



(60) 



The effective interaction strength F is assumed to be 
small: 



rz/< 1 



(61) 



The condition, eq.(|61|), allows to calculate correc- 
tions due to electron-electron interactions perturbatively. 
From calculations done at equilibrium it is known that 
while this model predicts the correct qualitative behav- 
ior for the electronic conductance, it misses the double 
logarithmic correction to the tunneling density of states. 
The different singularities in the conductance and the 
tunneling-density-of-states have been attributed to the 
presence (absence) of an approximate gauge invariance 
symmetrycj in the former (latter). Gauge invariant quan- 
tities are less sensitive to the details of the electron- 
electron interaction. Since the current correlation func- 
tion is evidently gauge invariant, we expect our model 
interaction to yield the correct (singular) temperature 
and voltage dependence, up to a non-universal numerical 
prcfactor. 

Following the rotation, eq. (p^ , of the interaction part 
of the action, eq.(59b), in Keldysh space, it reads 



iS['l>] = iFTr{$^(r,t)cri$(r,t)} . 



(62) 



Next one performs a gradient expansion of eq.(59c) 
treating the field (j) under the logarithm as a small per- 
turbation, which yields 
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iSl'P, Q] = -^Tr{I?(VQ)2-4i(0„7" + e)g} . (63) 

The correlation function of the fluctuations of the electric 
potential, corresponding to eq.(|6^), is 

{Mr, oj)Mr', -^)) = -irS{r - r')al'^^ , (64) 

and the correlation function of w and w are given by 
eq.(|3^). Let us now use this model to calculate the in- 
teraction correction to the d.c. value of electron current. 
Employing eq.(|3^) we obtain 

51 = ^ — X 

4 

(Tr{(QVQ- VQQ)72}Tr{0„7"g}Tr{0;37''Q})o , (65) 

where the subscript means that we perform averaging 
employing the correlators of eqs.(|6^ and (|35|). 

The analogous correction to the current noise is given 

by 

(0) = - ^M.,M.,,]. 



Tr{</.„7"Q}Tr{0^7'5Q} 



(66) 



So far our analysis was general and applied to diffusive 
conductors of any dimensionality. From this point on we 



are going to focus on quasi-two-dimensional systems and 
study eqs.(|65|) and (mb for such geometries. Evaluating 
the expression of eq.(|65|) (see Appendix we obtain a 
correction to the electron current. The strongest func- 
tional dependence of that correction on temperature (or 
voltage) is obtained when the parameter C,—T[i&yi{eV^ kT) 
is larger than the Thouless energy, i.e.a Q ^ Eth ■ Our 
analysis then yields 



^Ohm 

2^ 



ln(CT 



(67) 



Eq.(|67|) for the effect of interaction on the current at 
d = 2 is valid to the leading logarithmic correction. It 
features both the voltage and the temperature in a sym- 
metric way, namely as cutoffs which govern the singular- 
ity. For high temperatures, T ^ eV, the singularity is 
determined by T, which is consistent with the celebrated 
Altshuler-Aronov correction to the conductance. 

We now turn to the calculation of the interaction cor- 
rection to the current noise eq. (|66|). 

For frequencies smaller than the Thouless frequency 
the current-current correlation function does not depend 
on the choice of cross-section. Moreover, we also asr 
sume that T is much larger than the Thouless energyE3. 
To calculate the correction to the current noise due to 
electron-electron interaction one expands eq.(|6^) in the 
fields w,w. After averaging over these fields (see ap- 
pendix ^) we finally obtain 



5S'=^(0): 



Gohn 



2Tln (^l-t-eFcoth 



Eq. (|6|) is the main result of this section. Consider- 
ing the asymptotic high temperature behavior of eq.(68), 
one notes that dS{0) is given to leading order by the ex- 
pression for Nyquist noise. Corrections to the noise in 
this limit are-xelated to the interaction corrections to the 
conductanceEj (the latter is related to the Nyquist noise 
through the FDT), namely 



2 ^^M? 



(69) 



Eq. (|69|) is related to the correction to the mean current, 
eq.(§^). 

In the large voltage limit, the shot component of the 
noise is dominant, but interaction correction to it are still 
determined by the temperature 



63" 



e\V\\ii 



t) 



(70) 



Similarly to previous high temperature limit, here too the 
interaction correction to the noise is negative (cf . . As 



we can see the results for the current (67|) and for the cur- 
rent noise ( |70| ) are no more simply related to each other. 

The result we have obtained can be qualitatively un- 
derstood in the following way. The charge transfer 



1^ -mW? 



(68) 



r 



through disordered junction is a stochastic process. For 
non-interacting electrons the transfer of the charged par- 
ticles is described bji-.a time sequence with a binomial 
distribution functioned. The mean current has to do with 
the mean charge transferred through a cross-section over 
a given time interval, while the second cummulants of the 
current noise is related to the dispersion of the distribu- 
tion function. Electron-electron interactions render the 
electron motion more correlated. This affects the disper- 
sion more strongly than it affects the average value of the 
distribution, implying that the ongoing random process 
is no longer binomial. The deviation from the binomial 
distribution, though, is small in the inverse dimensionless 
conductance. 

Finally we would like to discuss in brief possible exper- 
imental implications of our results. The effect predicted 
here can, in principle, be tested experimentally in sam- 
ples where the conventional Altshuler-Aronov suppres- 
sion of the conductance is observed. The non-equilibrium 
character of the quantities considered raises the ques- 
tion whether other mechanisms (such as inelastic elec- 
tronic collisions) which might modify the noise spectrum 
as well, would not mask our effect. We thus would like to 
point out a qualitative difference between the interaction 
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corrections to the noise discussed here, and correctimis 
that arise because of the inelastic length being finitecJ'EJ. 
Normally the electron inelastic length increases as the 
temperature and/or the applied bias decrease. For this 
reason corrections associated with the ratio L/lin being 
finite diminish as one lowers the bias (temperature). In- 
deed, when the inelastic length is determined by large 
momentum transfer the correctipOptf) the shot noise due 
to inelastic collisions is given byulll 



el Utt^ K3 [eVLf 
y 40320 Dep 



0.02 • e^I 



(71) 



Here K3 = ^/lAire^v) is the inverse screening length in 
three dimensions. The second equality in eq.(f7l|) makes 
use of the expression for the three-dimensional inelastic 
length (large momentum transfer), given by 



Ij, 



64 



1 



{eVf 



(72) 



This expression is to be used for quasi-two-dimensional 
conductors (metallic films). While the leading term, eq. 
(p3|), is linear in the voltage, the corrections are propor- 
tional to its second power, and therefore can be neglected 
for small enough voltages. The corrections we have 
found, eqs.( |6^ , [70|) , behave in the opposite way. They 
become increasingly more pronounced as one lowers the 
voltage and the temperature. For a quantitative estimate 
we consider a metallic film of length L = 10~*cm and of 
thickness 10~^cm, having sheet resistance of 10 f2/square. 
The diffusion coefficient is taken to be Z? = lO'^cm^sec"^, 
the elastic mean free time - t = 10~^'^sec and the Fermi 
energy ep = 0-leV. We consider a voltage bias of 
V = lOfiV. The current noise in this range of param- 
eters is of the order of 10~^^Coulomb Ampere. The cor- 
rections to the current noise due to inelastic collisions 



are in this case of order 10~^^Coulomb Ampere, while 
the corrections due to the quantum interactions arc of 
order 10~^^Coulomb Ampere. 
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APPENDIX A: NOISE OF NON-INTERACTING 
ELECTRONS 

In this Appendix we derive eq.(^). One may note 
that there are two types of terms in eq.(^5|). The "dia- 
magnetic" part (arising from the term in the action) 
is proportional to the dimensionless conductance, and 
therefore it is sufficient to consider its value at the saddle 
point: 



[l-F{x,e-n)F{x,e)] 



dedfl exp(iili) 



(Al) 



The (MM) correlator is proportional to the square of the 
dimensionless conductance. It contains both reducible 
and irreducible contributions to the current-current cor- 
relation. In order to extract the contribution of this cor- 
relator to the second order cummulant of the current (ir- 
reducible contribution), we will need to expand it up to 
second os^r in w, w. Expanding M to first order in these 
variableala one obtains 



Mi{x,t) = J-^ j exp{int)dedn 



Vw(x, e, en)+Vuj{x, e, e - O) - V<Z>(x, e, e - n)F{x, e)F{x, e - Q) + 



\/F{x, e - n)F{x, e)w(x, e,e-Q) + VF{x, e)F{x, e-Q)u;{x, e, en) 



(2) 



( Here and later on the subscript k in the 1^ and M^ denotes the power of fiuctuating field up to which it is expanded). 
Now one needs to average the product of Afi's over the fluctuating fields 

(Tr{(QVQ - VQQ) 72}.,tTr{(QVQ - VQQ) 72}.',*')-,^ = 
2 



VV'Dfj(l - F,^,,F,_n,x') + VV'D_o(l - F,,,F,+o,..) + 



TTi^D J dxiVViDaV'ViD_n(l-i^...iF,_o.xJ 



(3) 



Combining eq.(^ and eq. (^) with eq.(p5[) one obtains the expression for the current correlation function through 
different cross sections at a finite frequency (for which the diffusion approximation holds) 
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2tt 



[1 - F{x, e)F{x, e - lu)]S{x - x') + 



S{x,x' ,uj) = e^TrvD J 



(4) 



The last result can be re-expressed in terms of an effec- 
tive electrons temperature (eq. (^)), leading to eq.(49) 



APPENDIX II: CORRECTIONS TO THE D.C. 
CURRENT 



In this appendix we derive eq.(67), i.e. the interac- 
tion corrections to the disorder-averaged value of the d.c. 
current. One first expands the operator M up to second 
J order in the fluctuating fields 



M2(a;i) - / d[e 



£3,^2 "I" '^a:i,ei,e3'^^a;i,e3,e2 



(1) 



where d[e\ implies integration over all e's, /de/27r. Next we evaluate the expansion up to second order in the fields 
w, w of the quantity 



(Tr{0i7iQ}Tr{0272Q})0 = / d[, 



Using the expressions (|65|,|l|;^D wc find that the correction to the d.c. current is given by 



(2) 



Wx 



(3) 



After averaging over the fields w and w we find that, to lowest order in the interaction amplitude and inverse 
conductance, the correction is given by 



61" 



-ieiT^V^DTA j d[e]dx2{Fx2,c3 - Fx2,e^)VFx^,eA^x^,X2,^3-<^^^X2,x^ 



(4) 



Now we will use the fact that the diffusion propagator at frequency f2 is a function of the difference of the spatial 
coordinates; it decays on the scale ^JD/Q, which is much smaller than the sample size, provided that TiVl ^ i?Th- This 
allows us to integrate over the relative coordinate (assuming that the distribution function does not change much on 
that distance), 



5P 



^e^VDTA^^^ f d[Se]d^[q]D^[q,-Se]coth^^^^^^ 



2ttL 



2kT 



(5) 



Integrating over q in (|^) we obtain 



f d6e , 


^Se + eV\ 


/-coth 




V 2kT J 



Restoring the "universal" ( i.e., F-independent) coef- 
ficient from the known equilibrium result we derive 
(6) eq.@. 



After integration over the energy Se we finally find the 
correction to the current which exhibits a logarithmic 
singularity. The latter is smeared on a scale which is de- 
termined by the maximal between the temperature and 
the voltage: 

e^TTiyTAV ^ ,^ , 



APPENDIX III: CORRECTIONS TO THE 
CURRENT NOISE 



In this Appendix we derive eq.(|68|). To find the in- 
teraction correction to the current noise we use the fact 
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that for frequencies lower than the Thouless energy the order respectively. We find that the contribution corn- 
result is independent of the cross-sections involved; it is ing from the corrections to MM vanishes after integra- 
therefore legitimate to integrate the current correlation tion over space coordinates (it contains terms of the type 



function over the volume of the system. 

Employing eq.(]6fl) we will need to expand the and 



^x,L — Dk,o)- To find the correction that arises from the 
"diamagnetic term" we need to expand the latter up to 



the MM terms in w and w up to a second and third second order in w and w, 
t') = f d[e] exp(z((6i - e^)t + (63 - £2)^')) 



2^ei,e2-fei-fe4^e3>e6^«6>e4 '^«3.e4''^ei ,£5 ,^2 "t" 2<5€3,e4 -fei -^£3 Wei ,65 "U^es ,£2 



(1) 



Combining eq.(|l]) with eq.(g) and using the fact that the diffusion propagator decays on scales much shorter than 
the system's size, we can convert the spatial integration over the relative coordinate to an integral in Fourier space ( 
analogously to how it was done in Appendix |l|). We obtain 



(5l£(c^) = ATA J d[e]d[q]dxlml^D^[q,ei - e2]|(^^x,ei - F,,,J(2^^,^,,F,,,,_„ - 1 



Thus the correction (eq. (|66|) to the noise is given by 



L2 



J d^[q]d[Se]lm^^D^[q,6e]^Y{Se,UJ,eV,T) 



where we have defined 



Integrating over energy and (" center-of-mass" ) coordinate we find that at zero frequency, lo — 0, 



Y{Se,uj,eV,T) = / d[u]dx{Fu,:, - F„_5.,x)(2F„,,F„_^,, - 1) 



(2) 



(3) 



(4) 



Y{Se,0,eV,T) 



2kT coth 



eV -Se 
2kT 



eV coth 



/ eV 



coth 



eV-5e 



coth 



coth 



eV + Se 



eV + Se 
2kT 



, , , ^„„.. , , , ,, - 2eycoth ( I coth f I 

y2kT J y V 2fcr J V 2fcT JJ \2kT J \2kT J 

Substituting eq.(||) into eq.(^) we conclude our calculation of the interaction correction to the noise 



5S'='=(0) 



GTiriy 
35(27r)3 



2fcrin(Cr) + el/coth 



/ eV 



\2kT 



(ln(Cr)+ln(fcTr)) 



(5) 



(6) 



Restoring a F- independent coefficient from the equilibrium calculation of (H), we finally derive eq.(68) 
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